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Abstract 

We present a new method to search for a light scalar top with < nit, decaying 
dominantly into a c-jet and the lightest neutralino, at the LHC. The principal idea 
is to exploit the Majorana nature of the gluino, leading to same-sign top quarks 
in events of gluino-pair production followed by gluino decays into top and stop. 

The resulting signature is 2 6-jets plus 2 same-sign leptons plus additional jets and 
missing energy. We perform a Monte Carlo simulation for a benchmark scenario, 
which is in agreement with the recent WMAP bound on the relic density of dark 
matter, and demonstrate that for nig < 900 GeV and nig > nig the signal can be 
extracted from the background. Moreover, we discuss the determination of the stop 
and gluino masses from the shape of invariant-mass distributions. The derivation 
of the shape formulae is also given. 

1 Introduction 

Owing to the large top Yukawa coupling, the supersymmetric partners of the top quark, 
the so-called ‘scalar tops’ or ‘stops’, play a special role in the MSSM [1]. This is manifest 
in i) the mixing of the left- and right-chiral states fL,R to mass eigenstates ti ^2 with 
the mixing and mass splitting proportional to Ml^ = 6.4^2 (A — ficot/3); ii) the large 
impact of the stop masses on the light Higgs mass mu through radiative corrections; and 
hi) the influence of the stop sector on the renormalization group (RG) running of the 
SUSY breaking parameters. RG running and L-R mixing can render the lighter stop, 
ti, much lighter than all other squarks. Indeed, there are scenarios which prefer a very 
light fi, lighter than the top quark. For example, the requirement of a strong enough 
hrst order phase transition to preserve the baryon asymmetry of the Universe, together 
with the Higgs mass bound from LEP2, favour < mt and somewhat small values of 
tan/3 ~ 2-8 [2-5]. A light ti also ameliorates the hue-tuning in SUSY models. 

The experimental bound on mi^ from LEP2 [6] is about 94-100 GeV, depending on the 
decay mode, the t mixing, and the mass difference between the ti and the lightest SUSY 
particle (LSP). This leaves us with 100 GeV < mi^ < mt as a very interesting possibility 
to consider. A ti in this range may be discovered at the Tevatron [7,8] provided the p- 
LSP mass difference is large enough. At the LHG, studies of scalar tops typically suffer 
from an enormous SUSY background, often from b decays, which makes it very difficult 
to extract the i signal [9-11]. However, it has to be noted that most of these studies 
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mi, [GeV] 

120 

130 

140 

150 

160 

170 

180 

Tevatron 

5.43 

3.44 

2.25 

1.50 

1.02 

0.71 

0.50 

a(tjj), LHG 

757 

532 

382 

280 

209 

158 

121 


Table 1: NLO cross sections in pb for ti pair production at the Tevatron and the LHC, 
computed with Prospino2 [12], For the radiative corrections, rriq = 1 TeV for all squarks 
apart from H, and rrig = 660 GeV. 

were performed within the mSUGRA scenario, which imposes relations between stop and 
sbottom masses, and also between stops/sbottoms and the rest of the spectrum. In this 
paper we follow a different approach, discussing the phenomenology of a light ti in the 
framework of the general MSSM. 

Quite generally, pair production of light stops has a large cross section at both the 
Tevatron and the LHG, comparable to the ti cross section. The cross sections at next- 
to-leading order (NLO) are given in Table 1. If < m^± -|- and mi^ — m^o < 
m\Y, the ii may decay dominantly into cy? [13]. This gives a signature of two c-jets 
plus missing transverse energy which may be extracted at the Tevatron provided 

ruji — m^o > 30 GeV and m^o < 100 GeV, for an integrated luminosity of 4 fb“^ [7,8]. 
At the LHG, however, it will be exceedingly difficult to use this signature for a discovery. 

In this paper, we therefore propose an alternative signature to search for a light ti at 
the LHG: For masses up to ~1 TeV, gluinos are also copiously pair-produced at the LHG. 
The cross sections at NLO are given in Table 2. If rrig > mi^ + rrit, gluinos decay into 
stops with a large branching ratio. The important point is that being Majorana particles, 
they decay into ttl or it^ combinations with equal probability. Pair-produced gluinos 
therefore give 

gg —> tii^il^ttilil, tti^i^ (1) 

and hence same-sign top quarks in half of the gluino-to-stop decays. Let now the W 
stemming from t bW in the same-sign top events decay leptonically, and let the ti 
decay into cXi- This gives a signature of two 6-jets plus two same-sign leptons plus jets 
plus missing transverse energy, 

gg bb 1+1+ (or bb ri-) + jets + , (2) 

which is quite peculiar. As we will show, it will serve to remove most of the backgrounds, 
both SM and SUSY, and may hence be used for discovery of a light stop at the LHG. It may 
moreover be used for determining a relationship between the gluino, stop and neutralino 
masses. Furthermore, it might offer a possibility to test the Majorana nature of the 
gluino by comparing the number same-sign and opposite-sign di-tops from Eq. (1). The 
signature Eq. (2) will therefore be useful for LHG analysis even in the case that the stop 
is discovered at the Tevatron. Of course, gluino-squark and squark-pair production with 
subsequent q ^ gq decay can also contribute to the signal, although with additional jets. 
For the masses of interest in this study, a{pp —> gq) is comparable in size to a{pp —> gg), 
while <j{pp —> qif*'^) is smaller by about a factor of hve. 

The rest of the paper is organised as follows. We hrst discuss in Section 2 constraints 
from the relic density of dark matter. In Section 3 we perform a case study of a light stop 
for the LHG: the parameters of our benchmark scenario are given in Section 3.1, the Monte 
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(^{gg) [pb] 

113 

31.6 

10.4 

3.84 

1.56 

0.68 

0.31 


Table 2; NLO cross sections in pb for gluino-pair production at the LHC, computed with 
Prospino2 [12]. For the radiative corrections, = 150 GeV and = 1 TeV has 

been assumed. 

Carlo simulation is explained in Section 3.2, the effective mass scale is shown in Section 3.3, 
the signal isolation is discussed in detail in Section 3.4, and the determination of masses 
in Section 3.5. In Section 4 we then present our conclusions. Finally, the derivation of 
the formulae for the invariant-mass distributions is given the Appendix. 

2 Constraints from relic density 

Requiring that the lightest SUSY particle (LSP) provide the right amount of cold dark 
matter 

0.0945 < < 0.1287 (3) 

at 2a [14,15] puts strong constraints on any SUSY scenario. In the standard approach, 
the relic density is oc 1/(av), where (av) is the thermally averaged cross section times 
the relative velocity of the LSP pair. This thermally averaged effective annihilation cross 
section includes a sum over all (co-)annihilation channels for the LSP. For a neutralino 
LSP, the value of hence depends on the Xi mass and decomposition (i.e. on the 
gaugino-higgsino mixing), as well as on the properties of all other sparticles that contribute 
to the annihilation and co-annihilation processes. The main channels are i) annihilation 
into fermion pairs via s-channel Z or Higgs exchange, ii) annihilation into fermion pairs 
via t-channel sfermion exchange, iii) annihilation into WW or ZZ via t-channel exchange 
of charginos or neutralinos, and iv) co-annihilation with sparticles which are close in mass 
to the LSP. The value of Qh'^ depends sensitively on the kinematics of the dominating 
process; in case i) e.g. on rriA — 2m^o, and in case iv) on the mass difference AM between 
the neutralino and the co-annihilating sparticle, often the lighter stau or as in our study 
also the light stop. As has been shown in [16,17], a shift in AM of only 1 GeV can induce 
an 0(10%) change in 

For a given set of gaugino-higgsino parameters one can hence derive constraints on 
[part of] the rest of the spectrum in order to satisfy the WMAP bound of Eq. (3). We 
illustrate this considering two scenarios motivated by the results on baryogenesis viable 
light stop models of [18], but neglecting GP-violating phases for simplicity. For computing 
the neutralino relic density, we use the program micrOMEGAs 1.3 [19]. In order not to 
vary too many parameters, we use a common mass scale of 250 GeV for all sleptons 
apart from fi, and a common mass of 1 TeV for all squarks apart from ti, assuming 
fi ~ tr and ii ~ i^. Varying and then means adjusting mp, m^. At, and A^ 
for hxed m?, u and tand. For the gaugino masses, we assume the GUT relation 
M2 = {g2/gi?M, ~ 2Mi. 








Figure 1: Neutralino relic density for Mi = 110 GeV, M 2 = 220 GeV, /r = 300 GeV, 
tan/3 = 7: in a) the WMAP allowed band in plane, and in b) VLh? as a function 

of rrti^ for various fi masses. Gomputed with micrOMEGAs 1.3. All masses are in [GeV]. 
The other parameters are as explained in the text. 

2.1 Ml = 110 GeV, = 300 GeV, tan/3 = 7 

As the hrst scenario we take the parameter point Mi = 110 GeV, M 2 = 220 GeV, 
H = 300 GeV and tan/3 = 7. This gives a spectrum of ^ = {105,191, 306, 340} GeV 
and = {189,340} GeV. The Xi is dominantly a bino with only 1% wino and 5% 
higgsino admixture. With slepton masses around 250 GeV, squark masses around 1 TeV 
and rriA S> 2m^o, the Xi annihilation cross section is much too low, leading to ~ 0.8, 
which is well above the WMAP bound. One possibility to achieve the right is to lower 
rriA to about 250 GeV. In this case the neutralinos annihilate efficiently through XiXi 
A —>■ bb, leading to Qh'^ ~ 0.1 (when closer to the pole of ruA — 2xi, the x*^’s annihilate 
too fast and Qh‘^ becomes too small). Another possibility is to rely on co-annihilation 
with stops or staus. This puts rather strong constraints on the fi and/or ii masses, since 
stau co-annihilation occurs for — m^o < 10 GeV while stop co-annihilation requires 
mi^ — m^o < 25 GeV. These constraints imply that the decay products of the staus or 
stops will be difficult to detect at a hadron collider. 

Figure la shows the WMAP-allowed band in the plane for ti ~ and 

Ti ~ tr, niA = 1 TeV and the other parameters as explained above. Figure lb shows the 
neutralino relic density as a function of mi^ for various values of mf^. As can be seen, 
the WMAP bound puts a lower limit on the ti and fi masses of mi^ > 125 GeV and 
ruf^ ^112 GeV (which is, however, only a stringent bound if one requires that the LSP 
provide all the cold dark matter). It can also put an upper limit on one mass as a function 
of the other. This bound is more severe. For example, for mi^ > 130 GeV agreement with 
WMAP requires ^114 GeV if no other process helps the Xi annihilate efficiently. 
The exact masses of the other sleptons and squarks are not important for the computation 
of the relic density in this example because they are too heavy to contribute signihcantly 
in (co-)annihilation processes. 
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Figure 2: WMAP allowed band in the plane analogous to Fig. la, but for p = 

180 GeV. 

2.2 Ml = 110 GeV, /x = 180 GeV, tan /3 = 7 

For the second scenario, we lower p to 180 GeV, keeping the other parameters as above. 
This leads to ^ = {93,152,188,274} GeV and = {139,273} GeV. The Xi has 
now 4% wino and 25% higgsino admixture and annihilates quite efficently into W~^W~, 
giving a relic density of = 0.138. We hence need only a small additional contribution 
to {(Jv), e.g. from light staus or stops, or from Higgs exchange. 

Figure 2 shows the WMAP allowed band in the plane for this scenario. The 

allowed region is much larger than in Fig. la, especially because for ^130 GeV, 
t-channel exchange of fi with nif^ < 160 GeV is sufhcient to bring into the desired 
range. For mf^ ^105 GeV, on the other hand, stau co-annihilation comes into play, driv¬ 
ing below the 2a WMAP bound. Likewise, turns out too low for my < 120 GeV 
because of co-annihilation with stops. Note, however, that in this scenario a ti heavier 
than about 144 GeV will decay into bxt rather than into cXi, which gives additional b 
jets and, more importantly, additional same-sign leptons. This changes the signature of 
Eq. (2) to 

gg bbbb /+/+/-/“ + , (4) 

or more generally to 45 -|- AW -|- It is clear that the information from the same-sign 

top quarks is lost in this way. 

3 Case study for the LHC 

3.1 Ghoice of parameters 

In order to explore the feasibility of extracting the light stop signal Eq. (2) at the LHG, 
we perform a case study for the parameters of Sect. 2.1 (Mi = 110 GeV, p = 300 GeV, 
tan/3 = 7), my = 150 GeV and rrig = 660 GeV, assuming BR(ti —> cXi) = 1- All 
squark mass parameters apart from my^ are set to 1 TeV in order to maximise BR(y — 
til) and suppress background from e.g. b —> tx~ decays. For the sleptons, we take 
My y = 250 GeV. As explained in Sect. 2.1, agreement with WMAP can be achieved 
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Table 3: Input parameters for the LSTl scenario [masses in GeV]. Unless stated otherwise, 
the SM masses are pole masses. 
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hi 
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660.00 

Xi 

104.81 

X2 

190.45 

Xs 

306.06 

X4 

340.80 

xf 

188.64 

X^ 

340.09 


h 

118.05 

H 

251.52 

A 

250.00 

262.45 




Table 4: SUSY mass spectrum [in GeV] for the LSTl scenario. For the squarks and 
sleptons, the hrst two generations have identical masses. 

by rriA — 250 GeV or, if is large, by ruf^ ~ 112-113 GeV. In what follows we use 
rriA = 250 GeV. We call this the LSTl scenario. The parameters and the resulting mass 
spectrum calculated with SuSpect 2.3 [21] are given in Tables 3 and 4, respectively. Note 
that the SUSY-breaking parameters in Table 3 are taken to be on-shell. 

3.2 Monte Carlo simulation 

We have generated SUSY events and tt background equivalent to 30 fb“^ of integrated 
luminosity with PYTHIA6.321 [22] and GTEQ 5L parton distribution functions [23], cor¬ 
responding to three years running of the LHG at low luminosity. The SUSY NLO cross 
sections are found in Table 5. We have also generated additional SM background in hve 
logarithmic pr bins from px = 50 GeV to 4000 GeV, consisting of 5 x 10^ of lU -|- jet, 
Z -I-jet, and WW/WZ/ZZ production and 3.5 x 10^ QGD 2 —> 2 events per bin. 

Detector simulation has been done using the generic LHG detector simulation 
AcerDET 1.0 [24]. This expresses identihcation and isolation of leptons and jets in terms 
of detector coordinates by azimuthal angle 0, pseudo-rapidity p and cone size AR = 
y^(A0p”+~(Ar^. We identify a lepton if px > 5(6) GeV and |? 7 | < 2.5 for electrons 
(muons). A lepton is isolated if it is at a distance AR > 0.4 from other leptons and jets, 
and if the transverse energy deposited in a cone AR = 0.2 around the lepton is less than 


































cr(titi) 



too 


a{qq) a{qq*) 

(^ixtg) 

a{tt) 

LSTl 

280 

5.39 

4.98 

1.48 

0.774 

0.666 0.281 

0.0894 

737 


Table 5: The main SUSY NLO cross sections in pb for the LSTl scenario, computed with 
Prospino2 [12]. For comparison, we also give the tt NLO cross section taken from [20]. 


10 GeV. Jets are reconstructed by a cone-based algorithm from clusters and are accepted 
if the jet has pr > 15 GeV in a cone Ai? = 0.4. The jets are re-calibrated using a flavour 
independent parametrisation, optimised to give a proper scale for the dijet decay of a 
light (100-120 GeV) Higgs boson. The 6-tagging efficiency and light jet rejection is set 
according to the pt parametrization for a low luminosity environment given in [25]. 


3.3 Effective mass 


In Fig. 3 we show the distribution of the effective mass Meff, 

+ (5) 

i 

for the LSTl scenario under study. The cuts used are: 

• Require at least four jets with > 100,100, 50, 50 GeV in each event. 

• No isolated electrons or muons. 


• > max(100,0.25^pJr ) GeV. 


The clear excess of events toward high values of Mgg, compared to the SM distribution, 
indicates that SUSY should be easily discovered in our light stop scenario. 

The effective mass has also been shown (see [26,27]) to work as a measurement of the 
effective SUSY mass scale, given by 


m?o 

= M -^ 

^'^^susy ^'^-'susy ^ 


susy 


( 6 ) 


where the SUSY mass scale, M^nsy, is a cross section weighted average of the masses of 
the initially produced SUSY particles; usually gluinos and squarks. In the light stop 
scenario, stop pair production overwhelmingly dominates the SUSY cross section. Yet it 
is difficult to extract the stop mass from the effective mass in this manner, since the cuts 
placed to handle the SM background also remove most of the stop pair production events. 
What remains is mostly gluino-pair production and gluino-squark production, and we can 
instead use the effective mass to estimate the gluino mass. In [27] the linear regression 
relationship Mgst = 1.7+ 134.15 GeV, where Mgst is the htted peak of the Mgg 
distribution, was found from a large set of random mSUGRA and MSSM models. The 
joint statistical and systematic error on the SUSY mass scale from using this relationship 
was estimated to be below 40% after one year running of the LHG at low luminosity. 
With the peak effective mass value of Mggt — 1023 GeV at LSTl, and a conservative 
estimate of accuracy at the level of 30% after 30 fb~^ of integrated luminosity, this gives 
^susy = 523 ± 157 GeV. This may serve as a hrst estimate of the gluino mass scale. 
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Figure 3: Effective mass distribution. SM contributions are ti (filled circles), PF+jet (tri¬ 
angles), Z-|-jet (inverted triangles), WWjWZjZZ production (stars) and QCD (squares). 
The sum of all SM events is shown by the hatched histogram. SUSY events are shown as 
open circles. 

3.4 Signal isolation 

The peculiar nature of the signal with a same-sign top pair will serve to remove most of 
the background, both SM and SUSY. We use the following cutsd 

• Require two same sign leptons (e or /i) with > 20 GeV. 

• Require at least four jets with > 50 GeV, at least two of which are 5-tagged. 

• > 100 GeV. 

• We cut on top content in the events by requiring two combinations of leptons and 
6 -jets to give invariant masses mu < 160 GeV, consistent with a top. 

The effects of these cuts are shown in Table 6. The column “21ep 4jet” gives the 
status after detector simulation and cuts on two reconstructed and isolated leptons and 
four reconstructed jets; “26” is the number of events left after the 6-jet cut, assuming a 
6 -tagging efficiency of 43%; is the cut on missing transverse energy and “SS” the 

requirement of two same-sign leptons. These cuts constitute the signature of Eq. (2). Note 
the central importance of the same-sign cut in removing the SM background, which at 
that point consists only of ti events. The cuts on transverse momentum and top content 
“2t” are used to further reduce the background. We hnd that the gluino-pair production, 
followed by gluino decay into top and stop and leptonic top decay, is easily separated from 
the background. 

^It should be noted that the cut values have not been optimized to isolate the signal for this particular 
model (LSTl); they have been chosen so as to effectively remove SM events and to isolate events with 
two same-sign top quarks. 
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jet 

Pt 
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77'miss 

2t 

SS 

Signal 









10839 

6317 

4158 

960 

806 

628 

330 

Background 








SUSY 

1406 

778 

236 

40 

33 

16 

5 

SM 

25.3M 

1.3M 

35977 

4809 

1787 

1653 

12 


Table 6: Number of events left after cumulative cuts for 30 fb ^ of integrated luminosity. 


To investigate other possible backgrounds to our signal we have used MadGraph II 
with the MadEvent event generator [28,29]. The search has been limited to parton level, 
as we hnd no processes that can contribute after placing appropriate cuts. We have 
looked at SM processes that can mimic a same-sign top pair by mis-tagging of jets or the 
production of one or more additional leptons, as well as inclusive production of same-sign 
top pairs. We assume that the two extra jets needed in some cases could be produced by 
ISR, FSR, or the underlying event. In particular we have looked at diffractive scattering 
qq — W^q'W^q' and the production of a top pair from gluon radiation in single W 
production qq' — tiW^. Also checked is the production of titt, titb, tibi, tW~tW~, 

iW^iW^ and W^W^Wj]. We place cuts on leptons and quarks as given above, and 
demand two lepton-quark pairs consistent with top decays. We also require neutrinos 
from the W decays to give the required missing energy. After these cuts and reasonable 
detector geometry cuts of AR > 0.4 and \ri\ < 2.5 for all leptons and quarks, we hnd that 
the cross-sections of all of these processes are too small, by at least an order of magnitude, 
to make a contribution at the integrated luminosity considered. 

Last but not least we have assumed that there is no additional same-sign top produc¬ 
tion from havour-changing neutral currents (FCNC), i.e. that the anomalous couplings in 
tgc(u) vertices are effectively zero. See [30] for a discussion on same-sign tops in FCNC 
scenarios. 

3.5 Determining masses 

Having isolated the decay chain it will be important to measure the properties of the 
sparticles involved to conhrm that the decay indeed involves a light scalar top. Since the 
neutralino and the neutrino in the top decay represent missing momentum and energy, 
reconstruction of a mass peak is impossible. The well studied alternative to this, see 
e.g. [26,31-34], is to use the invariant-mass distributions of the SM decay products. Their 
endpoints can be given in terms of the SUSY masses involved, and these equations can 
then in principle be solved to give the masses. 

In our scenario there are two main difficulties with this. First, there are four possible 
endpoints: and of which the hrst simply gives a relationship 

between the masses of the W and the top, and the second and third are linearly dependent, 
so that we are left with three unknown masses and only two equations. Second, because 
of the information lost with the escaping neutrino the distributions of interest all fall 
very gradually to zero. Determining exact endpoints in the presence of background, while 
taking into account smearing from the detector, effects of particle widths etc. will be very 
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Figure 4: Invariant-mass distributions for LSTl. These distributions only take into ac¬ 
count the kinematics of the decay, i.e. no spin or width effects are included. 


difficult. The shapes of the invariant-mass distributions are shown, for some arbitrary 
normalization, in Fig. 4. 

In what follows we have partially solved the second problem by extending the endpoint 
method and deriving analytic expressions for the shape of the invariant-mass distributions 
rribc and rriic- The derivation is given in Appendix A. These shapes can now be used to 
£t the whole distributions of the isolated signal events and not just the endpoints. This 
greatly reduces the uncertainty involved in endpoint determination, and provides the 
possibility of getting more information on the masses. One could also imagine extending 
this method to include spin effects in the distribution, to get a handle on the spins of the 
SUSY particles involved and possibly conhrming the scalar nature of the stop^. 

In htting the rribc and rriic distributions, we start from the isolated gg signal of Sec¬ 
tion 3.4. However, these events contain some where one or both of the W decay to a tau, 
which in turn decays leptonically. These taus are an additional, irreducible background 
to our distributions. The 6-jets and leptons are paired through the cut on two t quark 
candidates. A comparison with Monte Carlo truth information from the event generation 
shows that this works well in picking the right pairs. The issue which remains is to iden¬ 
tify the c-quark initiated jets and to assign these to the correct 6-jet and lepton pair. The 
precision of our mass determination is limited by systematics from these problems. 

Two different strategies can be used for picking the c-jets. The strong correlation 
between the tagging of 6- and c-jets suggest an inclusive 6/c-jet tagging of at least four 
jets per event. The two types of jets can then be separated on their 6-tagging likelihoods, 
and the requirement of two top candidates in the event. A thorough investigation of this 
strategy requires a full simulation study, using realistic 6-tagging routines. The other 
strategy, which we follow here, is to accept a low 6-tagging efficiency to pick two 6-jets 
and reject most c-jets. The likelihoods in the 6-tagging routine can then help pick the 
correct c-jets from the remaining jets. In this fast simulation study we are restricted to 
a simple statistical model of the efficiency of making this identihcation. We have looked 
at two cases: One worst case scenario with no direct c-jet identihcation, where we only 
use the kinematics of the event to pick out the c-jets, and one where we have assumed an 

^For more details on the derivation of invariant-mass distributions in cascade decays, and the inclusion 
of spin effects, see [35]. 
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additional 20% probability of identifying a c-jet directly from the 6-tagging likelihood. 

For events where we have missed one or both c-jets, they are picked as the two hardest 
remaining jets with < 100 GeV. The upper bound on transverse momentum is applied 
because the stop is expected to be relatively light if our signal exists, and it avoids 
picking jets from the decay of heavy squarks. Our c-jet candidates are paired to the top 
candidates by their angular separation in the lab frame, and by requiring consistency 
with the endpoints of the two invariant-mass distributions we are not looking at. E.g. 
if we wish to construct the mtc distribution we demand consistency with the endpoints 
and Events with no consistent combinations of c-jets and top candidates 

are rejected. 

The ht functions for and rriic are given in Eqs. (24) and (48). In principle both 
of the two linearly independent parameters and a could be determined by hts. 

However, we typically have mtmp m? for light stops, so that a ~ 1. For LSTl, the 
nominal value is a = 0.991. The distributions are sensitive to such values of a only at 
very low invariant masses. Because of the low number of events, no sensible value can be 
determined from a £t; we therefore set a = 1. The fit quality and value of is found 
to be insensitive to the choice of a for a > 0.980. 

r\j 

The results of the hts to assuming no c-jet tagging, are shown in Fig. 5. The 

combined result of the two distributions is = 383.2 ±4.9 GeV, to be compared with 
the nominal value of 391.1 GeV. The large values of the hts and the low value of 

indicate that there are some signihcant systematical errors. Comparing to Monte Carlo 
truth information we have found that the peaks at around 100 —150 GeV, responsible for 
the bad ht quality, are chiehy the result of events with one or more tans (see above). In 
Fig. 6 we show the results assuming 20% c-tagging efficiency. The combined result has 
improved to = 389.8 ± 5.3 GeV. We expect to be able to do better than this with 
full information from the 6-tagging routine. 


4 Conclusions 

We have investigated a baryogenesis-motivated scenario of a light stop (m^^ < m^), with 
ti cXi as the dominant decay mode. In this scenario, pair production of ti leads to 
a signature of two jets and missing transverse energy, which is all but promising for the 
discovery of ti at the LHC. 

We have hence proposed a method which uses instead stops stemming from gluino 
decays: in gluino pair production, the Majorana nature of the gluino leads to a peculiar 
signature of same-sign top quarks in half of the gluino-to-stop decays. For the case that 
all other squarks are heavier than the gluino, we have shown that the resulting signature 
of 26’s ± 2 same-sign leptons ± jets can be extracted from the background and 

serve as a discovery channel for a light ti. 

We have also demonstrated the measurement of a relationship between the gluino, stop 
and LSP masses. Taken together with a determination of other invariant-mass endpoints, 
and a measurement of the SUSY mass scale from the ehective mass scale of events, this 
may be sufficient to determine the masses of the SUSY particles involved. In particular, 

^We require that the values are below the rough estimates = 430 GeV, = 480 GeV and 

= 505 GeV, approximately 40 GeV above the nominal values, so no precise pre-determination of 
endpoints is assumed. 



Events/10 GeV/30 fb Events/10 GeV/30 fb 




Figure 5: Invariant-mass distributions without c-jet tagging (black with error bars) and 
best fit. Left panel shows mtc, right panel mic. Also shown are the contributions from 
the SM background (green) and the SUSY background (blue). The SUSY background 
consists mostly of events with one or more tans (see text). 




Figure 6: Invariant-mass distributions with 20% c-tagging efficiency after 6-tagging. For 
details see caption of Figure 5. 


if the invariant-mass distributions of the isolated events £t the predicted shapes, this 
strengthens the interpretation of the events as gluino decays into a top and a stop. 

Last but not least some comments are in order on the robustness of the signal. While 
our analysis has been done for rrig = 660 GeV, we have checked that the signal remains 
significant enough for a 5a discovery for gluino masses up to rrig ~ 900 GeV. In the 
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high-mass region, gluino-squark production with subsequent squark decay to gluino is 
responsible for a significant fraction of events. On the other hand, if sbottoms are lighter 
than assumed for LSTl, they will contribute to the SUSY background through hi —>• tx^ 
{i,j = 1,2) decays. For instance, lowering the sbottom masses by a factor of 2, while 
keeping all other parameters as in LSTl, would increase the background after cuts by 
5 events. While this is not dramatic, the fits to invariant-mass distributions are worse 
because of the opening up of new gluino decays g —>■ hhi. Lowering the stop mass while 
keeping the LSP mass fixed, entering the stop-coannihilation region, should also make 
the signal more difficult to find because the c-jets become too soft. Still, we find only a 
12% decrease in the signal when setting mi^ = 120 GeV. This can be explained by a large 
number of signal events passing the cuts because of squark decays to gluinos producing 
an extra jet. This implies that the same-sign tops can be used to search for a light stop 
even in the stop-coannihilation channel, unreachable at the Tevatron. The fit to the 
invariant-mass shapes is again worse than for LSTl. 

We conclude that if ti decays into cy?, light stops may be discovered through the search 
for same-sign tops in the decay of pairs of gluinos for a wide range of SUSY masses. A 
relation between the gluino, stop and LSP masses can be determined from invariant-mass 
distributions. In this paper we have demonstrated the feasibility of our method; a full 
simulation will clearly be necessary to assess its full potential. 
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A Derivation of shape formulae 


The differential decay width for two independent angular variables (cosines) u and n, in 
a decay with a given spin configuration, is 

f = fn{u)e{u - 1)0(1 - u)f^{v)e{v - 1)0(1 - v), (7) 


where the spin information is contained in the functions /«(«) and fv{v) and the theta 
functions are ordinary step function, limiting the values of u and v that give non-zero 
decay width. 

By a change of variables one can easily go from an expression of the invariant-mass in 
terms of these variables, to the distribution sought. Given the invariant-mass m{u, v), we 
have 


1 aP 

P dm'^ 



1 d^T 
fdm^dv^^ 
d{u,v) 1 d'^T 
d (m^, v) P dudv 


dv. 


An extension to three or more independent angular variables is trivial. 


(8) 









We now proceed to derive the shape of the distribution of the invariant-masse of the b- 
and the c-quark, and the invariant-mass of the lepton and the c-quark, mic- We will 
assume that all particles are spin-0, so that we have isotropic decays in particle rest frames. 
This lets us set /„ and to be constant. This result is what we expect when summing 
over all hnal states. Effects of non-zero spin can be introduced via these functions. We 
also assume that the lighter quarks, b and c, are massless. 


A.l rribc 


In the rest frame (RF) of the stop, energy and momentum conservation in the decay of 
the gluino, stop and top gives 


Eh 





2 2 

2 2 2 
m~ — mf — m- 

y ^ t\ 

2mi^ 

_ ml - rri^ _ 

2El' - 2^(E‘'Y - mj cos 0%' 


(9) 

( 10 ) 

( 11 ) 


where E^}, and E^^ are the energies of c, t and b respectively in the stop RF and 6\l 
is the angle between t and b in that frame. 

We can rewrite cosOh in terms of an isotropically distributed angle, in the top 
RF.^ By expressing in the two rest frames, and since cos 9^1 = cos9^-\, we get 


d‘(d‘ - V(4)" - = E‘(Et - ^(E'0^ - m|cos«!,). (12) 

In the top RF, from the conservation of energy and momentum in the decays of the gluino 
and the top, we have that: 


^ ml - 
= 


El = 


2mt 

m? -|- ml — ml^ 

2mt 


(13) 

(14) 


Solving for cosOh, using (10), (11), (13) and (14), we then arrive at 


a — cos 9i 


gb 


cos 6**1 = 

1 — a cos 9i 


gb 


where a is given by 




a = 


m 


2 ’ 
1 


(15) 


(16) 


^We here make the assumption that the top has spin-0. Again this is what a sum over final states will 
yield. 












with mi and m 2 defined as 


2 2 2 2 
mi = rUg — , 


4 A A 2 2 

m 2 = mi — 4m^ m^^. 

We can now find an expression for : 


(17) 

(18) 


m 


be 


= 2El^Eb i 1 


= [m 


cosell 


0 , 1 — a cos 6 *h 

max\4 I go 


be 


1. Qj 


cosell 


(19) 


where the endpoint of the distribution, is 

, {ml - ml,) H - ™|;) (™? + ”>2) 


(m 


be 


mt 


Using the switch of variables 


u = 


1 - cos 


V = 


2m~ 

Cl 


1 -COS0ffe 


we can write the invariant-mass as 


ml = imZn' 


1 — a -|- 2av 

1 -|- fl 


u. 


( 20 ) 


( 21 ) 


( 22 ) 


From Eq. ( 8 ), and setting fu{u) = fv{v) = 1, we find the distribution of the invariant-mass 

e (u) e {1 — u)e (v) e {i — v) dv 

1 „ /. ml{l + a) 


1 <9F 

poo 

d {u, v) 

F dml 

J —00 

d {ml, v) 


1 a 


{m^rrJo 1 


a 2av 


e 1 


(1 - a + 2av) 


dv. (23) 


The theta function leads to two different lower limits on the integration, depending on 
the value of mbe- Performing the integration gives 


1 -|- (3 


1 or 

fdml 


1 -|- (3 


2a{mfl^^)^ ■■■ 1 - a 


In 


‘'be 

1 -|- (3 


for 0 < ml < 


2a mi^ 


^max'|2 (24) 

In-^^—^— for ^ < ml < 


expressing the distribution in terms of two parameters, a and m^^^. 

A .2 mic 

In the top RF, energy and momentum conservation in the decay of the gluino, stop, top 
and W gives, in addition to Eq. (13), 


El = 


mtm^ 


{ml + ml) - {ml - m^) cos e{^i ’ 
mt{m^ 


(25) 


■0 ™i;) 


m\ — m 2 cos 0)^^ ’ 


(26) 






















where Ej and E^. are the energies of b, I and c respectively in the top RF, 9\^i is the angle 
between W and / in that frame and is the angle between ti and c. 

We again change angles to isotropically distributed angles (under spin-0 assumptions). 
From cos 9 \yi = — cos6*^;, and expressing mh in both the top and W rest frames, we have 

ElEHl + cosely,) = £“'<(1 - cosO. (27) 


In the W RF, energy and momentum conservation in the top and W decays gives 


eV = 


EY = -TTlw- 


m: 


m 


w 


2mw 


(28) 

(29) 


Using (13), (26), (28) and (29) we hnd that 


cos 9 \yi — 


{ml - m\y) - {ml + m^) cos 9^ 
{ml + mly) - {ml - m^) cos 9^' 


(30) 


From cos = cos 6*|^, and expressing in both the t and ti rest frames we get 


gc 


E'M - = Et(4 - \/(d‘)'“ ™i"1) 

from which, using (9), (10), (14) and (26), it follows that 


(31) 


, a - cos 

cos9j= -^ 

^ 1 — a cos 6*h 


(32) 


gc 


This could also have been found from (15), using a symmetry argument. 
We can now hnd an expression for ml^ : 

ml = 2EjEl{l-cos9l) 

1 — a cos \ f ml + mly — {ml — m^) cos 9^ 


= 


where 


1 “h Cl 


{rn^y = 


2ml 


l-cos0f. 


(33) 


{mi + Tnl) 


Using the switch of variables 


u = 


1 -cos0^ 


V = 


2m? 

ll 


1 - coses' 


(34) 


1 — cos 9i 


W = 


Ic 


(35) 


we can write 


ml = (m—) 


max 42 / 1 - a + 2aM\ /m^ + (m^ - m^) v 


1 -j- n 


mt 


w. 


(36) 
















Shape 

y{xi) 

y{x2) 

I 

y2 < y{xi) 

2/1 < yix2) < 2/2 

II 

yi < y{xi) < 2/2 

2/1 < yix2) < 2/2 

III 

1/2 < y{xi) 

y{x2) < yi 

IV 

2/1 < y{xi) < 2/2 

y{x2) < 2/1 

V 

y{xi) < 2/1 

y{x2) < 2/1 


Table 7: Shapes for the area of integration for the integral in Eq. (41). 


Since 






r dudvdmf^ 
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(1 + a)9{u)6{v)6{w) 


(^max )2 ("x _ a + 2au) (m^ + {mj - 
where 6{x) = 6{x)6{l — x), we have from a generalization of (8) that 
1 9r m^(l + a) 6{w) 
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mw)v] 


vY 


dudv. 


(37) 


(38) 


Jn Jo (1 - (I + 2a'(()(m^ + (m) 

The step functions for w give more complicated bounds on the integration than what was 
the case for mbc- We must have > 0 and 

m‘f{l + a)ml 


m 


w 


+ (m; - m^)v > 


(^max)2(l_Q + 2an)' 


(39) 


This integration and the integration limits from the last inequality are best explored by 
yet another change of variables 


X = 1 — a + 2au, y = + (m^ 

from which the integral can be written 
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1 ar 
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1 “h Cl 
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Ic 
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' 1—a 


[m 


Ic 
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dxdy. 


(40) 


(41) 


We note that there is a maximum of rrif^ < as expected. There are then Eve 

different possible shapes for the areas of integration in the xy—plane. For ease of notation 
we use 

xi = l- a, X2 = l + a, yi = y2 = m^. (42) 

The shapes can be categorized by bounds on y{xi) and y{x 2 ) as found in Table 7. The 
inequalities in the categorization of the shapes can be expressed as bounds on mf^: 


yi 
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(43) 

yi 

< 

y{x2) ^ (m--)2^ < ml 
mf 

(44) 

y2 
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y{xi) ^ (m;);“)2 ^ J < ml 

(45) 

y2 

< 

y{x2) ^ {rnH'^Y < ml 

(46) 






















We can see that there will be two cases, depending on the mass hierarchy: 


Case A: 


Case B: 


1^ 




< 
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1—a 

1+a 

1—a 

l-\-a 


( 47 ) 


We give the invariant-mass distribution for both of these cases. 


Case A Here Shape III is excluded, so the invariant-mass distribution is 
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(48) 


Case B Here Shape H is excluded, giving 
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